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Abstract We report measurements of coherence times of an rf SQUID qubit using pulsed 
microwaves and rapid flux pulses. The modified rf SQUID, described by an double-well 
potential, has independent, in situ, controls for the tilt and barrier height of the potential. 
The decay of coherent oscillations is dominated by the lifetime of the excited state and low 
frequency flux noise and is consistent with independent measurement of these quantities 
obtained by microwave spectroscopy, resonant tunneling between fluxoid wells and decay 
of the excited state. The oscillation's waveform is compared to analytical results obtained 
for finite decay rates and detuning and averaged over low frequency flux noise. 
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1 Introduction 

In the last few years there has been significant progress towards improved coherence times 
in superconducting qubits [1,2,3,4]. However, progress has been slower in rf SQUID qubits 
then other similar devices such as phase qubits [2][3) and persistent current qubits [5]|4]|6). 
This lack of improved coherence times is usually attributed to the larger size of rf SQUID 
devices which require a large loop (i.e. 150 fim x 150 /xm) to provide the necessary ge- 
ometrical inductance. This paper will report recent measurements of coherence times of 
a modified rf SQUID, including Rabi oscillations, that suggest that the dominant source 
of decoherence in these devices is low frequency flux noise from two level fluctuators as 
seen in other superconducting qubits. We propose taking advantage of the sensitivity of rf 
SQUIDs to low frequency flux noise to investigate the source of the flux noise in the local 
environment. 

An rf SQUID qubit in its simplest form consists of a superconducting ring of inductance 
L interrupted by a single Josephson junction with critical current l c shunted by a capacitance 
C(see Fig.[T]inset). When a static external flux (& x ) is applied to the SQUID loop, it induces 
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Fig. 1 (Color online) The potential of an rf-SQUID at a j3 = 1.32 and V = 0.505 showing localized energy 
levels and the corresponding value of mean flux (green x). The inset shows a schematic representation of an 
rf SQUID 

a screening current /, = —I c sin(2^4"/4>o) where <t>o = ^ is the flux quantum. The screened 
flux linking the loop (<P) must satisfy, <P = <P X +LI S . The dynamics of this device is analo- 
gous to a particle of mass C with kinetic energy jC<t> 2 moving in a potential given by the 
sum of the loop's magnetic energy and the Josephson junctions coupling energy. Expressing 
the fluxes in units of <Pq, this potential is (7J 

C/(<?>,W = C/o[2^ 2 ^-^) 2 -jS L cos(2^)] , (1) 

where Uq = t^/Ai^L, Pi = 2tiLI c /4>q. Lowercase (<f>__) is used to denote the various fluxes 
in units of <Pq (e.g. (j) x = <P x /<Pq). For 1 < j3i < 4.6 and <f> x near 1 /2, the potential, shown in 
Fig.[U has the form of a double well with the two wells representing different fluxoid states 
of the rf SQUID and having screening currents (of equal magnitude at the symmetry point, 
tf) x = 1/2) circulating in opposite directions around the loop [17 j. Changing (j) x has the effect 
of tilting the double well potential. One can select either the lowest level in each well or the 
two lowest levels in a single well for the basis states of the qubit. Since the two states have an 
easily measurable flux difference in the double well basis, this is referred to as the flux basis 
and requires a relatively low barrier to provide adequate coupling between the states. The 
single well mode, where the basis states are coupled by microwave radiation, on the other 
hand, requires a relatively high barrier so that the basis states do not couple to the states in 
the other well. This mode of operation is very close the that of qubits commonly referred 
to as "phase" qubits, so we will use this term to identify it. Figure [2] shows an example of 
calculated energy levels for the potential as a function of ty x . As can be seen, changes in (j> x 
cause a much greater change (about 100 x) in the level spacing for the flux basis than for 
the phase basis. Since, as we shall see, there is substantial noise in (j) x , this gives a strong 
advantage to the phase basis in terms of coherence times. 
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Fig. 2 (Color online) The energy levels of the rf SQUID referenced to the lowest energy level in the left well. 
The blue dashed line (magenta dotted line) is the first (second) excited state localized in the left well. The 
black lines are the levels localized in the right well. 

2 Design and Fabrication of the Qubit and On-Chip Control and Read-Out Circuitry 

The qubits used in this work, shown schematically in the bottom center of Fig. [3ja), are 
enhanced versions of the simple qubit discussed above. This design uses a gradiometric 
configuration to help isolate the qubit from fluctuations in the ambient flux on the length 
scale (150 fim) of the qubit, but can still be described by the simple potential of Eq. [T] 
If a single flux quantum were trapped in the outer loop when the qubit cooled through its 
transition temperature, then the motion of the system between the two wells of the potential 
would correspond to motion of the flux from the upper to lower loop of the gradiometer 
through the junctions connecting these two loops. If the ambient flux were zero, this situation 
would also result in the qubit potential being symmetric without further bias flux. 

As noted above, different modes of operation of the qubit requires changes in the barrier 
height. It would, therefore, be very convenient to be able to adjust the height of the barrier 
of the double well potential, in situ, independently of (f> x . This requires the ability to vary 
the critical current of the Josephson junction coupling the two wells. The twin requirements 
of low capacitance and rapid modulation using small control currents mean that a single 
large area junction cannot be used for this purpose. However, the critical current can be 
adjusted if the single junction is replaced by two junctions, with critical currents I c \ and 
I c 2, in a small loop, a dc-SQUID, with inductance I as shown in Fig. [3]a). By applying a 
flux, fyxdc, to this small loop the critical current of the dc-SQUID is suppressed allowing 
the barrier of the double well potential to be modulated. The addition of a second junction 
formally makes the potential two dimensional. However, if L 2> £ and l c \ w 7 e 2, then this 
2D potential is very well approximated by the one dimensional potential in Eq.Q]wifh small 
corrections |8 9| which can be neglected. jSl is now given by Pl{^>xcIc) — j8locos(^^ c ) 
where jS L0 = 2nL(l c \ +I c2 )/$o. 

The flux (j) x that tilts the potential has both a low frequency component, <j) x £f, that is 
applied using the on-chip gradient coil shown at the bottom right of Fig. [S^a) and high 
frequency components, discussed below in sec. III. The gradiometric design of the qubit has 
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Fig. 3 (Color online) (a) Schematic and (b) micrograph of an rf SQUID qubit and the readout magnetometer, 
(c) a cross section of the wafer around the junctions and (d) micrograph giving a detailed view of junctions. 



the further advantage that, to first order, there is no cross-coupling between the tilting flux, 
<f>xlf, and the barrier modulation flux, ty x d c - The outer loop of the 150 jUm x 150 /xm qubit is 
a 5 \im wide Nb film with a thickness of 150 nm as is the film separating the two halves of 
the gradiometer except the thickness is 200 nm. This gives an effective qubit inductance of 
L m 215 pH. The dc-SQUID loop is 5 jUm x 5 |im giving £ w 10 pH. These inductances are 
calculated using the 3D-MLSI software package [ 10] and are consistent with measurements. 

The fluxoid state of the qubit is measured using a hysteretic dc SQUID magnetometer 
shown at the top of Fig.[3ja). The measured mutual inductance between the qubit and the 
magnetometer is 5 pH, consistent with the design value of 4.3 pH. The self inductance of the 
magnetometer is 56 pH and the junctions are 2.15 x 2.15 jtm and 2.85 x 2.85 /im giving 
a maximum critical current of 10.8 fiA. This asymmetric design, along with the asymmetry 
in the magnetometer flux bias loop, fa, is part of a scheme to decouple the magnetometer 
bias leads from the qubit to first order while permitting fast triggering of the readout. 

These circuits are fabricated at Stony Brook in our laboratory using a niobium trilayer 
process called the"SAL-EBL" process 111 II 12 1 . This process has been developed from our 
well established planarized process for Nb superconducting electronics - the PARTS-EBL 
process 1131 . The SAL-EBL process was developed specifically to study effects of various 
process steps on coherence in Nb based qubits. In particular, to maximize junction and 
film qualitv ll 1 II 141 and to minimize the turn-around time only essential process steps were 
retained, eliminating CMP (chemical mechanical polishing) and reducing process steps such 
as reactive ion etches (RIE) known to reduce film quality under certain conditions 1 15']. The 
resulting process has only one RIE step and uses a self-aligned lift-off (SAL) of the dielectric 
to define the junctions. By using lift-off for majority of the patterning, the process allows 
material flexibility for dielectric and wiring metal, and by using a combination of deep UV 
and electron beam lithography (EBL) on the same resist layer at each step allows great 
design flexibility. 
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The devices are fabricated on 2" oxidized Si wafers. The Nb/A10 x /Nb trilayer is pat- 
terned via lift-off. Both Nb and Al are dc sputtered in Ar and A10 x is formed by thermal 
oxidation. The junctions defined by EBL on a negative resist are formed by RIE in SFg gas 
with a minimum device size of 0.15 x 0.15 Mm 2 . The high area uniformity of the junctions 
allow us to accurately control junction asymmetry. The same resist mask is used to lift off 
rf sputtered SiC>2 dielectric. A Nb wiring layer is also patterned through lift-off. Fig. [3jc) 
shows a cross-section of a completed device indicating various layer thicknesses. The pro- 
cess turn-around time is approximately one week, and the use of EBL for all critical features 
enables rapid design changes. 

The size for the qubit junctions vary depending on desired qubit parameters and is 1.45 x 
1.45 jUm 2 for the data shown. The designed current density for these devices is 1 jixA/jUm 2 . 
Both the junction size and current density are measured on diagnostic chips and agree with 
designed values. In particular the measurements of electrical size from both current and 
resistance scaling are consistent with physical size of the junctions estimated from RIE 
undercut. The asymmetry of the junctions (1.6%) is extracted from fits to (j> as function of 
<j> x and is small enough to justify the use of Eq.|T][ 16]. 



3 Apparatus and Measurement Techniques 

Two essential criteria for the measurement of qubit dynamics are that one be able to prepare 
and readout states of the qubit with a time resolution short compared to the decoherence 
time of the states being studied and that the qubit be very well isolated from its environment 
so that the decoherence time is not affected, e.g. by the measurement and control processes, 
while the qubit is undergoing coherent evolution. This section presents the approaches used 
to achieve these goals with nanosecond resolution for qubits using the phase basis, which is 
far less sensitive to low frequency noise in (j) x . 

The on-chip bias coils, discussed in Sec. [2] can apply fluxes of the order of <Pq but are 
limited to frequencies of less than 100 kHz. To achieve nanosecond timing resolution, a sep- 
arate set of high frequency flux bias signals, a microwave ((j) xr f) and a video ((j) xp ) pulse with 
lengths on the order of nanoseconds and sub-nanosecond rise times, are inductively coupled 
to the qubit via a hole in the ground plane of a well characterized microstrip transmission 
line [ TTj on a separate chip suspended 25 fim above the qubit chip . These signals, while 
very fast, are limited in amplitude to several thousandth of <Pq. This gives, for the total bias 
flux applied to the large loop of the qubit, (j) x = (j) x (f + <j) xr f + (j) xp + (j> ff, where (j)„ff is a 
constant that includes the ambient flux and the zero offset in (j) x e/, which is taken to be 
at the symmetry point of the potential. The desired state of the qubit is obtained from the 
initial ground state of the left well by applying one or more microwave and video pulses. 

Since the magnetometer cannot distinguish between the two basis states in the phase 
basis, an intermediate step must be used for the readout of the qubit state. To accomplish 
this, at the moment readout is desired, a short pulse, (j> xp , is applied to tilt the potential 
rapidly but adiabatically to reduce the tunneling barrier between wells. The barrier height 
during this tilt is about 1.1 K for the excited state of the qubit giving a tunneling rate to 
the right well of about 10 9 s _1 while the rate for the ground state is about a factor of 500 
slower. Thus it is possible to clearly distinguish these two states by appropriate selection 
of the pulse length. After the end of the pulse, the system will be trapped in either the left 
well (ground state) or the right well (excited state) for at least 1 s. Subsequent readout of the 
occupied well using the magnetometer will then reveal the state of the qubit when <j) xp was 
applied. 



Pulse n 
Generator _l L 




^ u-strip ground plane 



Fig. 4 (Color online) Circuit schematic including the various stages of filtering and the components used to 
create and combine tj> u -f and <j> x „. 



The components that create and combine the microwave pulses and readout pulse are 
shown at the top of Fig. [4] The pulse generator is capable of producing measurement pulses 
with rise times as short as 200 ps. Two microwave mixers are used to modulate the envelope 
of the continuous microwave signal produced by the microwave source giving an on/off ratio 
of 10 3 . The output of the mixers is then amplified 20 dB and coupled through variable at- 
tenuators that set the amplitude of the microwave pulses applied to the qubit. Finally, these 
microwave pulses are combined with the video pulses using a hybrid coupler. The video 
pulses enter on the directly coupled port while the microwave pulse are coupled using the 
indirectly coupled port. The video pulses used for the data shown in the following sections 
are around typically 5 ns with a rise time of 0.5 ns. The directionality and frequency re- 
sponse of the coupler allows the two signals to be combined with a minimum amount of 
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reflections and loss of power. This signal is then coupled to the qubit through a coaxial line 
that is filtered by a series of attenuators at 1.4 K (20 dB), 600 mK (10 dB) and at the qubit 
temperature (10 dB) followed by a lossy microstrip filter that cuts off around 1 GHz. 

The low frequency bias and readout circuits for the qubit and magnetometer also have 
been optimized to control and measure the qubit while protecting it from noise. All signal 
paths are designed to be symmetric with respect to the chip to reduce the effect of common 
mode noise, and grounded signal sources are coupled through battery powered unity gain 
isolation amplifiers that effectively disconnect these signals from the earth ground. These 
low frequency lines are then filtered using EMI filters on entering the cryostat so that the 
dewar itself acts as an rf-shield for the cold portion of the experiment. They are further 
filtered using cascaded, 4-stage RC filters anchored at 1.4 K. Finally the lines are filtered 
using specially designed lossy microstrip filters at the qubit temperature. 

Common RC filters are limited to about 1 GHz due to the parasitic inductance of the 
capacitor. The lossy microstrip filters were designed to cut off exponentially below 1 GHz 
and are effective to much higher frequencies. They consist of a thin film of chromium de- 
posited on a glass or sapphire substrate and diced into long narrow chips. When the chips 
are placed in a brass housing, which acts as a ground plane, they form a lossy microstrip 
transmission line. The filters that must pass dc current also have a Nb meander line to pass 
the low frequency control currents without heating the filter assembly. These filters should 
be effective to the lowest order waveguide mode of the brass housing at 63 GHz while all 
other modes are greater then 500 GHz. At 18 GHz the microstrip filters have a measured 
attenuation greater than 95 dB. 

The qubit chip and all the necessary wires and cables are housed in a superconducting 
NbTi sample cell that acts as a magnetic shield below its transition temperature around 10 
K. This sample cell is mounted in the bottom section of a vacuum tight filter can, which has 
two chambers separated by an if tight block housing the microstrip filters. This filter can 
is placed inside a double layer Cryoperm magnetic shield on the sample stage of a dilution 
refrigerator capable of reaching a loaded base temperature of 5 mK. 



4 Results and Analysis 

There are a number of different ways of characterizing the coherence times in our flux qubit. 
Experimentally, the simplest method of probing the decoherence involves measuring the 
lineshape of macroscopic resonant tunneling (MRT) between fluxoid states [18|. However 
the dynamics that determine the tunneling process are complicated, especially when the 
width of the peaks are not much smaller than the spacing between them. The measurement 
flux pulse discussed in the previous section makes it possible to measure the occupancy of 
the excited state and study intrawell dynamics. This along with controlled microwave pulses 
makes possible a direct measurement of the lifetime of the excited state and the observation 
of Rabi oscillations. 



4.1 Lifetime of the Excited State 

The decay rate between the first excited state and the ground state in the same well (F) 
provides an upper limit on the coherence time of the qubit in the phase basis. F can be mea- 
sured directly by measuring the occupation of the excited state, pn, as a function of time. 
For this measurement, the qubit is pumped to the first excited state using a long microwave 
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pulse followed by a measurement pulse (j)_ xp to read out the occupation of this state as a 
function of the delay, % delay , between the end of the microwave pulse and the measurement 
pulse. A 80 ns microwave pulse, which is much longer than Tl = 1/Fi , is used in order to 
insure a constant initial mixture of the ground and excited state. 

Figure [5] shows the occupation of the excited state as a function of delay between the 
end of the microwave pulse and the readout pulse. The microwave frequency of 17.9 GHz 
is resonant with the energy difference between the ground and excited state for the value 
Qxtf = 0.050 x 10~ 3 used. The solid line is the fit to the data for pn °c exp(— Fjf) giving 
F = 0.05 ns~' or Tl = 20 ns. As we shall see below, such good fits are not always obtained 
for arbitrary bias. 

The lifetime of the excited state is similar to lifetimes measured in a number of other 
Nb/A10 x /Nb qubits operated in the phase basis (energy eigenstates in the same potential 
well). Paik et al. have observed decay times in dc SQUID phase qubits [ 19 1 around 13 ns for 
a Nb/A10 x /Nb junctions and 20 ns for A1/A10 X /A1 junctions. Martinis et al. observed a 
lifetime of 20 ns in a large area current biased Nb/A10 x /Nb Josephson junction [20]. They 
were able to improve the lifetime to 500 ns [3 1 through various changes including switching 
to Al junctions, reducing the size of the junction (~ 10 /im 2 ), and changing to SiN x for the 
wiring insulation layer instead of Si02- 



4.2 Level Spectroscopy 

If the decay with lifetime T 1 were the dominant source of decoherence in the qubit, one 
would expect a resonance in the microwave occupation of the first excited state given, in the 
low power limit, by p\\ « 2 r y/(<5 2 + Fj 2 /4) where 8 is the mean difference between the 
photon energy and the level spacing. For the potential parameters used for the resonant data 
shown in Fig.|6]the dependence of 8 on 4> v ;y is 8 = 1.70 x 10 9 s _1 /m4>o- Measurement of 
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the equilibrium occupation of the excited state near resonance is made using a readout pulse 
immediately following a 100 ns microwave pulse, which is much longer than the lifetime 
of the excited state as measured in Fig. |5j ensuring that the qubit reaches an equilibrium 
mixture of the ground and excited state. The detuning around the resonance is accomplished 
by changing <j>xef, rather than the frequency. This eliminates the effects of slight variations 
in (j) xr t as a function of the frequency. Figure [6] shows the occupation of the excited state 
for this potential as a function of 8 for 3 different microwave power levels spaced by 3db, 
thus giving a factor of two total variation in <j) xr f. For these values, the linewidth is nearly 
independent of <f) xr f. So, one would expect the data to be in the low power limit with a 
linewidth given by F . As can be seen, however, the actual linewidth is about an order of 
magnitude greater, indicating additional noise broadening. 

The effects of extra noise on this resonance between the ground and excited states within 
one well as well as on the resonant flux tunneling rate between levels in different wells (see 
Sec. |4.3l > can both be described by a similar set of evolution equations for the density matrix. 

In describing both types of resonant transitions between two levels, one needs to sepa- 
rate two types of interactions with the environmental noise. One is the decay, with the rate 
F, of the excited state to the lower-energy states within the same well induced by the noise 
components at frequencies ~ co p . Another is due to the noise components at frequencies be- 
low (Op, which induce fluctuations of the energy difference v(f) between the resonant levels. 
In the case of tunneling between the opposite wells, these fluctuations are induced directly 
by the energy shifts of one well relative to the other. For transitions within one well, both 
levels have practically the same value of average flux and the noise-induced fluctuations of 
the energy difference v are much smaller, roughly by two orders of magnitude for the param- 
eters of our qubits, than for the states in the opposite wells. We assume that the temperature 
T of the environment is sufficiently large (this assumption is made more precise later in this 
Section, and in Sec. 14.4b , so that the fluctuations v(f) induced by the low-frequency part of 
the noise can be treated classically. In this case, evolution equations for the elements py of 
the density matrix in the basis of the two resonant states that account for both relaxation and 
the low-frequency noise can be written as (see, e.g., [21 ]) 

Pn =^(Pio-Poi)-Fpn, 

Poo = '^(Poi-Pio)+Fpn, (2) 
Poi = -iv(t)poi - !^(poo — Pit) - yPoi , 

pio = ''v(f)pio + ^(poo-Pii)- yPio- 

The coupling amplitude a of the two states stands here for the tunnel amplitude A in the case 
of flux tunneling between the two wells, or for the microwave amplitude, (j) xr f, expressed as 
the Rabi frequency Q on resonance in the weak relaxation/dephasing limit (Sec. 14.41 ). 

For microwave excitation of the first excited state in the well, F is equal to the relaxation 
rate F of this first excited state, and v(f) in Eqs. (O is the detuning between the excitation 
energy of this state and microwave frequency. If the low-frequency noise is negligible, v(t) 
reduces to the average detuning 8, and the steady-state occupation of the excited states is 
obtained from the stationary solution of the Eqs. (f2]l and has a Lorentzian line shape as a 
function of 8 : ^ 

Pn{5) = 48i + mi + rf (3) 
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If the cut-off frequency CO c of the flux noise is much lower than the actual linewidth W of 
the resonant transition, determined self-consistently by the rms value u of the magnitude of 
the low-frequency noise and by the relaxation rate F (see Eq. © below), the noise can be 
accounted for simply by convoluting the intrinsic Lorentzian lineshape ® with some static 
distribution of the detuning values. Under a natural assumption of the Gaussian noise, this 
gives 

<pn> = -i- (dve-^l^puiv). (4) 

\/ZKU J 



In principle, the shape of the resonant peaks can provide information not only on the 
magnitude of noise but also on its spectral properties. For this, one needs to solve Eqs. (O for 
the time-dependent v(f) = 8 + v(f). This can be done analytically for the small microwave 
amplitude, Q. <C W, which can be treated as perturbation, i.e., when the contribution of Q to 
the peak broadening <[3j> is negligible. Spectral properties of the stationary classical Gaussian 
noise are characterized by the correlator: 

(v(f 1 )v(f 2 )) = -L f d(0S((Q)e i( °("-^ . (5) 



2k 

Perturbation theory in the coupling amplitude can be done directly in Eqs. $2% by integrating 
the equation for one of the off-diagonal elements of p under the assumption that diagonal 
elements are constant, substituting then the result in the equations for the diagonal elements, 
and finally averaging the resulting transition rate over the Gaussian fluctuations of detuning 
using Eq. ©. This straightforward calculation gives 

X2 2 f°° 

(Pn) = TF ^Te" r ' T / 2 cos5Tx 



2Ti Jo 

r 1 f , sin 2 cot/2, 
ex P{ / dfflS(ffl) 5-^-}. (6) 

7t J (0 A 

This result is the limit of the classical noise, reached for large temperatures, T 2> W, of a 
more general quantum expression obtained in |22|. Qualitatively, the quantum component 
of the noise leads to a finite shift of the resonant peak position away from the zero-detuning 
point 5 = [23 1. The shift disappears with increasing temperature T . 
For noise with a Lorentzian spectrum: 



<Pn) = ^=r { 
21 1 Jo 



dxe r ' T//2 cos5T 



x 



Eq. I© gives 



This equation describes the transition from the regime of the broadband noise, co c ^> W, 
when the noise simply increases the linewidth W of the Lorentzian (without the Q in the 
denominator) to Fi + 2u 2 /co c , to the regime of the quasistatic noise (O c <C W, when the 
lineshape is given by Eq. (O, and for strong noise, u ^> F, is: 

(Pii) = ^^=-^ 6/2 " • (9) 
^ l \2ku 



11 



1 1 1 1 1 1 1 1 1 1 1 1 

-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 

Detuning (Angular Frequency 10 9 ) 

Fig. 6 (Color online) The occupation of the excited state as a function of detuning for microwave powers 
corresponding to attenuator settings of 39 (squares), 36 (circles), 33 dB (triangles). Lines are fits using Eq. 
|4]at microwave amplitudes corresponding to the measured Rabi frequency for each attenuator setting (0.017, 
0.024, 0.034 ns ) with F = 0.055 ns~' convoluted with static Gaussian noise with o"V = 0.235 ns~' at the 
angular' frequencies of the Rabi oscillations that correspond to these microwave powers 

In the case of frequently assumed 1// noise of amplitude A: 

5(a)) = fi (10) 

Eq. ^ always gives the quasistatic result Q (with logarithmic accuracy) with u 2 ~ [AA/n) ln(W /(D/), 
where CO/ is the lower cutoffs of the spectrum d 1 0b roughly given by the time of one mea- 
surement cycle. 

The lines in Fig. [6] are fits to the data for the microwave-induced occupation of the ex- 
cited state, using Eq.[4]assuming extra low frequency flux noise noise in 8 of rms amplitude 
u = <7 V to obtain the lineshape. The best fit to the data, for all the powers shown, occurs for 
CT V = 0.235 ±0.01 ns~' and F = 0.055 ±0.005 ns~'. This analysis is based on the assump- 
tion that <j) x remains constant during the approximately 100 ns time of a single measurement 
(the quasi-static approximation) suggesting that low frequency flux noise is indeed the dom- 
inant source of resonant broadening. As a test of the degree to which this fit actually tells 
us something about the frequency dependence of the noise, we refit the lowest peak in Fig. 
fusing Eq.[8] which assumes Lorentzian noise with and arbitrary cutoff, co c , but is valid 
only in the low power limit. The results of this are that the data can be fit almost as well 
for values of co c up to 2 x 10 9 s _1 by increasing rms noise to <7 V = 0.28 ±0.01 ns~'. So it 
would seem that the fit provides very little information on the frequency dependence of the 
noise but gives a rather firm estimate of its amplitude. 

These fits can be extended to higher rf amplitudes where power broadening is important 
using Eq. [4] as shown in Fig. 13a), which shows the maximum occupation of the excited 
state as a function of microwave amplitude for a much wider range of (j) xr f. The red solid 
line is a fit to the peak amplitudes (5 = 0) using the same fitting parameters as in Fig. 
[6] For comparison, the blue dotted line is calculated for no flux noise, i.e. a v = 0, and 
F = 0.055 ns while the purple dashed line is calculated for <7 V = and F = 0.20 ns~' 
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Fig. 7 (Color online) (a) The occupation of the excited state on resonance versus microwave amplitude in 
units of the corresponding Rabi frequency. The lines are calculations for the following parameters; red solid 
F = 0.055 ns~' and o"V = 0.235 ns -1 , blue dotted line o~ v = with F = 0.055 ns -1 , purple dashed line 
Og — ar >d r = 0.20 ns . (b) The width of the spectroscopic peak from the Gaussian fits as a function of 
microwave amplitude in units of Rabi frequency. The lines are calculations for the following parameters; red 
solid F = 0.055 ns~> and ff v = 0.235 ns 1 , green dashed line o" v = with F = 0.055 ns 1 , blue dotted line 
o- v = and T = 0.585 ns-'. 



which gives the best fit to the low power data without considering low frequency noise. 
The peak amplitude cannot be fit over the whole range of microwave amplitudes without 
including the low frequency noise. Figure [T^b) shows the width of the peaks as a function 
of (j) xr f for the same data and the red solid line the fit for the same parameters as above. The 
green dashed line is the width for F = 0.055 ns -1 and the blue dotted line for F = 0.585 
ns~' and no low frequency noise. Again, low frequency Gaussian noise is required in order 
to obtain a reasonable fit to the overall lineshape of the data. 

As noted above, the data shown in Figs. [6] and [7] were obtained by changing <j) x (f over 
a narrow range near points where the spectrum was well behaved. In general, however, the 
spectrum is much more complex as can be seen in Fig. [8] showing the occupation of the 
excited state as function of both frequency and <j) x if for j3i = 1.30 for a much wider range 
of parameters. The color contours are proportional to occupation of the excited state, blue 
being lower occupation and red being the highest occupation (0.325). The curving ridge 
running diagonally from the top left to the bottom right corresponds the resonance between 
the ground and excited state. The black solid line shows the calculated difference between 
these states as a function of <j) x , where L and C are fitting parameters consistent with simula- 
tions and measurements (not shown) of thermal escape rates, positions of resonant tunneling 
peaks and photon assisted tunneling peaks 1161 . The gaps in the excited state occupation 
along this line are consistent with coupling to two level fluctuators, e.g. in the dielectrics, 
that have also been observed by other groups investigating superconducting phase qubits 
[24,3,251- In general these gaps do not correspond to the bias values where the states in the 
left well align with excited states in the right (indicated by the black circles). 

Another prominent feature in these data is the strong resonance near 17.25 GHz over 
a wide range of <j>xif, which we attribute to a cavity resonance in the NbTi cell. The larger 
peak of this cavity resonance, to the left of the resonance between the ground and excited 
state, corresponds to a two photon resonance to the second excited state of the qubit. The 
frequency of this cavity resonance is in the range of the lowest order cavity modes of the 
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Fig. 8 (Color online) The measured occupation of the excited state after a long microwave pulse expressed as 
color contours (blue being lowest and red being the highest) as function of both frequency and V . The solid 
lines are calculations of the energy level spitting between consecutive eigenstates states in the same well for 
P = 1.30, L = 190 pH and C = 210 fF. 
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Fig. 9 a) The measured occupation of the excited state after a long microwave pulse expressed as color 
contours (blue being lowest and red being the highest) as function <j) x and the delay between the end of the 
microwave pulse and the beginning of the readout pulse, b) Slices of the contour plot at (f) x (f = 0.767 x 10~ 3 



(black squares), 2.078 x 10 
the eye. 



(blue triangles) and 5.464 x 10 (red circles). The lines are just to help guide 



sample cell, even though most of these modes should be suppressed due to the extremely 
small dimensions (~ 1 mm) of the cell out of the plane of the chips. 

These anomalous behaviors are further illustrated in Fig. HJa), which shows the occu- 
pation of the excited state, using color contours, as a function of tdeiay and tyxtf- Here, the 
microwave frequency has been adjusted for each value of <j) x (f in order to remain on reso- 
nance. The bias value used for the measurement of the lifetime of the excited state shown in 
Fig.[5]is at the extreme left side of this plot. 
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At a number of bias points, the decay of the excited states is non-exponential, in contrast 
to the data of Fig. [5] For example, in the region near the gaps in the excitation spectrum, 
there are actually peaks in the occupation of the excited state as a function of t delay Figure 
gib) shows slices of Fig.^a) at faf = 0.767 x 1(T 3 (black squares), 2.078 x 10~ 3 (blue 
triangles) and 5.464 x 10~ 3 (red circles). At <j) x ef = 0.767 x 10~ 3 the decay is similar to 
that of Fig. [5] while for the other two <j> x values the behavior is clearly not exponential. The 
occupation of the excited state actually increases with time and reaches a local maximum at 
a non zero delay. At these bias points, it appears that the qubit is coupled strongly to a two 
level fiuctuator as was reported by Cooper et al. 1261 in their investigation on current-biased 
phase qubits. They also observed an oscillation superimposed on the decay of the excited 
state near one of these splittings in the excitation spectrum. 



4.3 Comparison with Resonant Tunneling Peaks 

As noted above the same equations (with different parameters) that determine the shape 
of the spectroscopic peaks discussed above should also describe the resonance tunneling 
between levels in different wells of the potential as initially seen in experiments on MRT, 
where an enhancement of the escape rate from one well to the other was observed when the 
lowest energy level of the initial well approximately aligns with an energy level in the other. 
When the energy difference between the levels (e) in opposite wells is much less than co p 
and the intrawell energy relaxation from the final state (the w th level of the right well) is r„, 
this tunneling rate, which has a form analogous to Eq.[3] is given by 1211 

r = (in 

Away from resonance when the levels are localized in their respective wells, the escape 
rate can be approximated as the sum of the tunneling rates from the ground state in the left 
well (|/)) to the lower energy states in the right well (|/)). Assuming Ohmic dampening, this 
gives for the individual rates 1271 

%, f = Y^R ( Ei ~ E A 1 ( '' /} |2 ' and Fesc = ^ 7 -f (12) 

where R is the ohmic dampening, Rq = h/4e 2 is the quantum unit of resistance and <t> is the 
flux operator. 

These rates can be measured by pulsing ty x from some initial value, where r esc in neg- 
ligible to the tilt at which r esc is to be measured and leaving it there for some time t m , in 
which the system has some reasonable probability of tunneling though the barrier into an 
excited state of the other flux well. Under the assumption above of incoherent tunneling, the 
system will then decay to the ground state of the new fluxoid well and remain trapped, so 
many repetitions of this measurement will give the probability of tunneling. The escape rate 
is determined from this probability of a transition, P(<j> x ), by 

r ^-h*(i=ks)- <I3) 

The measured escape rates are shown in Fig. [lO]as a function of fat. The resonant 
tunneling peaks are labeled with an index that indicates which level in the right well is 
aligned with the ground state in the left well. The red solid line is a calculation of the escape 
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Fig. 10 (Color online) Measured escape rates from the ground state in the left well as a function of V at a 
(S = 1.412. The solid red line is a calculation for Fj2 = 0.68 ns - '. The green dash dot line corresponds to the 
solid line with the addition of static Gaussian noise with o~ v = 0.23 ns -1 . The blue dashed line has the same 
level of Gaussian noise but with F12 = 10.5ns -1 . 



rates using Eqs.QT] and Q~2] with the measured qubit parameters. The decay rate, F12 = 0.68 
ns -1 is the product of the decay rate directly measured in Section |4~T1 and the level number 
(n = 12). This calculation gives peaks much narrower than the measured peaks. The green 
dash dot lines are these same peaks convoluted with static Gaussian flux noise in analogy 
to Eq. |4]with a v = 0.235 ± 0.01 ns , the same as the level of flux noise used to fit the 
spectroscopy in Sec. 14.21 For these parameters, the shape of the resonant tunneling peaks 
are well fit, but not the intervening valleys. In order to make both the peaks and valleys 
fit, the relaxation rate must be increased to F12 = 10.5 ns -1 which is inconsistent with the 
measured Fj . It may well be that a more refined theory is needed for the tunneling minima 
between the peaks. 



4.4 Time domain measurements 

While the data in the previous sections, clearly show the effects of various processes leading 
to decoherence in the qubit and even provide a means to obtain qualitative estimates of the 
decoherence rates, they are all obtained from non-coherent processes. In order to demon- 
strate coherence in our qubit and to compare the resulting decoherence rates with those due 
to the know sources of decoherence discussed above, we now turn to time domain measure- 
ments of the evolution of coherent superpositions of the two qubit states. 

An analytic solution to the Bloch equations for the population of the excited state in 
the presence of CW microwave, assuming no decoherence and that the system starts in the 
ground state at time zero, is given by 

q2 

P11M) = ^T^§2 ^ 2 (Q Rabi (8)t/2), (14) 
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where 



£i Rabi {8) = Vs 2 +n 2 . 



(15) 



This solution, shows that when driven with microwaves, the population of the excited 
states oscillates in time, demonstrating the phenomenon of Rabi oscillations, which have 
been seen in a number of different superconducting qubits [28 4 29 20 30]. Q, the fre- 
quency of the oscillations for 8 = 0, is ideally proportional to the amplitude of the mi- 
crowaves excitation <j> xr f . This linear dependence is accurately seen in our qubit at low 
power levels (see Fig. QT]upper inset) providing a convenient means to calibrate the ampli- 
tude of the microwaves incident on the qubit in terms of Q . 

The measurement sequence for coherent oscillations is very similar to that used to mea- 
sure the lifetime of the excited state except that the duration of the (generally shorter) 
microwave pulse is varied and the microwave pulse is immediately followed by the read- 
out flux pulse. This signal, on the high frequency line, is illustrated in the lower inset in 
Figure [TT| Figure QT| shows an example of the Rabi oscillations when 5 = and the mi- 
crowave frequency, f xr f = 17.9GHz. This bias point lies in a "clean" range of the spectrum 
(17.6 — 18GHz) as seen in Fig. [8] which should be the best region for observing coherent 
oscillations between the ground and excited states. Most of the time domain data, including 
the lifetime measurements of Fig. [5] have been taken in this frequency range. Each data 
point in Fig. [8] corresponds to the average of several thousand measurements for a given 
pulse length. 

4.5 Rabi Oscillations 

In order to simulate the time dependent dynamics of the Rabi oscillations in the presence 
of the low-frequency noise, it is convenient to separate this noise in the two parts. One, 
the quasi-static part, consists on the noise components with frequencies smaller than the 
linewidth W of the corresponding resonance. Under the conditions of the Rabi oscilla- 
tions, W is smaller than the strength Q of coherent rf coupling between the two levels, 
which gives the Rabi oscillation frequency in the case of resonance (8 = 0) and weak re- 
laxation/dephasing. This part of the noise can be taken into account by simple averaging 
over the static distribution of detuning 8 as in Eq. |4] The second, "dynamic" part of the 
low-frequency noise has components with frequencies extending beyond W. We treat this 
noise making an assumption (relevant for the data presented in this work) that the tempera- 
ture T is larger than not only W, but also the Rabi frequency Q. In this case, it is possible 
to characterize this part of the noise as the classical noise with the same constant spectral 
density 5 extending above Q and equal to the spectral density at frequency Q . As discussed 
in Sec. 14.21 such a broadband noise simply increase the decoherence rate of the two states 
by 5/2 (i.e., by u 2 /a c in the example of Eq. dHJ of that Section) to y = (T +5)/2. With this 
modification of the decoherence rate, Eqs. (|2]l for the evolution of the density matrix of the 
two rf-coupled states can be written as 



in terms of the imbalance p = (pn — poo)/2 of the occupation probabilities poo,Pn of 
the ground and excited state, and the real/imaginary parts m, v of the off-diagonal element 



P 



u 



Qv-Tp-r/2, 

Sv — yu , 

—Su — Qp — yv; 



(16) 
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Poi = u + iv of the density matrix. As Bloch equations dl6l > describe the dynamics of the 
Rabi oscillations, and are appropriate when the relaxation/dephasing rates are small on the 
scale of plasma frequency. The finite relaxation/dephasing rates make the analytic solution 
to these equations, which is shown below, considerably more complex than the undamped 
solution of Eq. 1141 

Decaying Rabi oscillations described by Eqs. l !16l > lead to some stationary values of 
the occupation imbalance p and coherences u, v that are established after the decay of the 
oscillations. The analytic solution of Eqs. dl6l > that corresponds to this process can be written 
as: 

A(f)=i4* + « -r XeVA;, (17) 

i 

where the "vectors" A are composed of the elements d 1 6b of the density matrix, A = {p, u, v}, 
and A st represents the stationary part of the solution: 

A ^ 8 2 + ff Q 2 7/r {S 2 + Y 2 ,ny,n5}. (18) 

The sum in Eq. H7[ is taken over the three eigenvalues Xj of the evolution matrix S: 

(y-T £2\ 
5= 5, (19) 




which are given by the roots of the cubic equation det(5— X) = 0. In the range of parameters 
relevant for the present discussion (not very strong relaxation/dephasing) one root is real: 

X l =X=s+-s- + (y-r)/3, (20) 

and the other two are complex conjugated: 

V3 

X 2 = A 3 * = K + ico, o) = — (21) 

K=5(j--s + ) + i(7-r), (22) 



where 



and 



s±= [(<? 3 +r 2 ) 1 / 2 ±r]' /3 , (23) 



i(i2 2 + 5 2 )-i(r- 7 ) 2 , 

r -r/i^ 2 S 2 : 



n z -8 z --(r-y) z j (24) 

The vectors Aj inEq. d 1 7b are the projections of the vector A(0) — A st of the deviations of 
the initial elements of the density matrix from their stationary values, onto the eigenvectors 
of the evolution matrix d 1 S>b that correspond to the eigenvalues Xj. These projections can be 
found conveniently from the following expression: 

*>= t-%l~-t m - A - h <25) 
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which for each j, directly cancels out components of the two other eigenvector k, I 7^ j in 
the initial vector A (0) — A st . 

If the dynamics of the system starts at t = in the ground state, then A (0) = { — 1/2,0,0}. 
In this case, using Eq. d 1 8b for the stationary vector A st and the matrix S d 1 9b in the Eq. | |25| |. 
together with the expressions for the eigenvalues Xj given above, we find the time depen- 
dence of the occupation probability pn of the excited state as 

7 Q. 2 

Pn - 



2F 8 2 + y 2 + Q 2 y/r 



(X - fc) 2 + co 2 



je A '(fc 2 + co 2 - 2fcy+ 7(7- F) - Q 2 - 8 2 r/y) + 

e a [(A (A - 2k) +2kj- y{y-r)+n 2 + 5 2 r/y) x 
cos cot+ [a> 2 (y-X) + (k-X)(kX -y(K + X)+ 

yiy-D-^-^r/y)^]} . 



(26) 



Equation l !26l > gives the decay rate of the Rabi oscillations as r\ = 7— fc. If the decay/relaxation 
rates in Eqs. fl!6b are small in comparison to Q, the rate r\ reduces to the Rabi oscillations 
decay rate F2 obtained previously by Ithier et al. |T), who expressed it in terms of the noise- 
induced decay rate at the Rabi frequency F v and the decoherence rate F, as 



ri 1 



8 2 +n 2 



£(3i2 2 + 25 2 )+5 2 F, + ^ 



(27) 



As discussed above, at large temperature T, the noise is classical at the Rabi frequency, and 
one can take F v = F^,. In this case, 7 = F/2 + F v , and for weak relaxation/dephasing, the 
two Rabi decay rates agree, 77 (5) = F2(5). 

The data for resonant Rabi oscillations, shown in Fig.[TTl The solid line is a fit using Eq. 
|26]including a 0.5 ns time delay for the rise time for the microwave pulse and averaged over 
quasi-static Gaussian noise in <j) x equal to that obtained from the fits to the spectroscopy data 
in Fig. [6] This gives a Rabi frequency /,„/,, = 119MHz and decay rate, 77 = 0.042 ns -1 , for 
the oscillations. Equation [27j with 5 = 0, together with our previously measured value for 
the decay rate of the first excited state F (which gives the rate F in the equations above), 
and the observed Rabi decay rate rj discussed below, imply F v = 0.01 ns -1 . Even though 
&Rabi is not affected by flux noise to first order for 5 = 0, the amplitude of the flux noise in 
our qubit is large enough to cause a measurable effect. If the flux noise were neglected, it 
would be necessary to increase Tj to 0.058 ns~' in order to account for the observed decay. 
In addition to increasing the decay rate, the low frequency noise also reduces the steady 
state occupation of the excited state, i.e. its value for long pulses, as calculated from Eq.|4] 
The ordinate of Fig.[6]has been calibrated using this value and is consistent with the estimate 
obtained from the calculated tunneling rates from the two levels involved during the readout 
pulse. 

The upper inset to Figure QT| shows the Rabi frequency extracted from Rabi oscillation 
measurements versus microwave amplitude in arbitrary units. As noted above, at low powers 
the Rabi frequency is linear in microwave amplitude as expected. At higher powers, how- 
ever, it begins to saturate. There is evidence that at these higher power levels there is a small 
probability of excitation to the second excited state. However since the second excited state 
has orders of magnitude higher escape rate, this small probability can be significant. 
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Fig. 11 (Color online) The occupation of the excited state as a function of the length of the microwave pulse 
demonstrating Rabi oscillations. The line is a fit to Eq. l26l for 5 = averaged over quasi-static noise with 
<y v = 0.22 ns~'. This gives f m ti = 119 MHz and decay time T2 = 24 ns~'. The upper inset shows Rabi 
frequency as a function of amplitude of applied microwaves in arbitrary units. The line is a linear fit to the 
lower microwave amplitude data. The lower inset show the measurement pulse sequence. 
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Fig. 12 (Color online) Rabi oscillations for detunings going from top to bottom of 0.094, 0.21 1, 0.328, 0.562, 
0.796 and 1.269 ns~' with the corresponding fits using F = 0.075 ns~' and 0" v = 0.22ns~' 



20 




2 4 6 8 10 12 2 4 6 8 10 12 

Time Between n/2 pulses (ns) Time Between nil Pulses (ns) 

Fig. 13 (Color online) (a) Measurement pulse sequence and (b) Measured occupation of the excited state 
after a series of 3 ns microwave pulses as a function of delay between the two pulses for various detuning 
pulse amplitudes (corresponding frequency); V HTs (no detuning), 0.1 V o's (0.21 ns -1 ), 0.45 V A's (0.96 
ns" 1 ), 0.8 V V's (1.70 ns" 1 ), 1.15 V o's (2.45 ns~') and 1.50 V <'s (3.20 ns" 1 ). (c) Calculations including 
noise for the pulse sequence by numerically solving the Bloch equations for the parameters corresponding to 
the measurement. 

The effects of flux noise on the decay of the Rabi oscillation, of course, becomes much 
more pronounced for 8^0. Figure [12] shows such Rabi oscillations at various detunings 
over a range of < 8 < 10 9 s~ l . At long microwave pulse times (not shown) the occupation 
of the excited state reaches the equilibrium values discussed in Sec l4.2l Near resonance, the 
value of pn for long pulse lengths is determined mainly by the noise amplitude, while for 
large 8 it is set by 8. The data cannot be fit to solutions of the Bloch equations using only the 
phenomenological decay constants. As in the spectroscopy, static Gaussian detuning noise 
must be included to fit the data over the whole range of detuning. The lines in Figure [T2l 
correspond to calculations with an initial Rabi frequency of 0.59 ns -1 , F = 0.075 ns~' and 
o"V = 0.22 ns -1 . The detuning frequency is taken from $ x using the conversion from Fig. 
[8] For the range of 8 in Fig. ll2l /l£ni is in a region with relatively few spectral anomalies. 
However for 8 < 0, on the other side of resonance towards splittings in the spectroscopy, 
the data generally do not agree with theory suggesting that the spectroscopic splittings have 
a strong effect on the coherence, as expected. 



4.6 Ramsey Pulse Sequence 

The effects of dephasing due low frequency flux noise should be seen most dramatically in 
the decay of Ramsey fringes. However the observation of Ramsey fringes in this system is 
complicated by the short lifetime of the excited state. To reduce these problems we have 
used the pulse sequence shown in Fig. 1131 a). Here the n/2 pulses are applied on resonance 
and the detuning is produced in between these pulses, using an additional (j> xp pulse to tilt 
the potential and slightly change the level spacing for a variable time. 

Figure \l3\ b) shows the measured occupation of the excited state as a function of delay 
between 3 ns microwave pulses (17.9 GHz) for several different values of detuning. The lines 
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are provided to help guide the eye between the points. The expected Ramsey fringes can be 
calculated by numerically solving the Bloch equations, using a time dependent detuning 
along with the time dependent driving amplitude. As in previous sections, the resulting 
oscillations are convoluted with quasi-static Gaussian distributed noise. Figure \l3\ c) shows 
these calculated occupations of the excited state for the parameters corresponding to the data 
in Fig. ll3t b). The scale factor for the amplitude of the detuning pulses was used as a fitting 
parameter but is consistent with the calibration of the detuning pulses from the ratio of the 
voltages between the measurement flux pulse and the detuning pulse. The rise times for the 
detuning pulse was assumed to be the same as was found for the measurement pulse. The 
rms value of the low frequency noise, a v , observed in Sec. l4.2l is used for this calculation, 
and the decay rate F that gave the best fit to the data is consistent with that measured in 
Sec. 14. II The determination of Fj has a relatively large uncertainty due to uncertainties in 
the exact shape and length of the tz/2 pulses, since the short coherence times are pushing 
the limits of our apparatus. However, the calculation shows qualitative agreement with the 
data. The limitations of the measurement do not allow for meaningful extraction of other 
parameters such as F^, . 

5 Sources of decoherence 

There is now a strong consensus that much of the decoherence observed in superconduct- 
ing qubits, including the flux noise seen here, comes from two level fluctuators (TLF) in 
dielectrics and oxides. The evidence seems to point towards something in the local environ- 
ment as the source of the decoherence. This is corroborated by observations of flux noise of 
unknown origin in SQUIDs 1 3 1 ,32,33] and persistent current qubits 1341 . 

The noise is typically independent of temperature T at low temperatures T < 0.1 K. 
This temperature independence of the noise, combined with the fact that the noise is pro- 
duced by an equilibrium source 1231 . indicates that the magnetic susceptibility of this source 
has paramagnetic l/T dependence on temperature [35,23]. While this observation supports 
the theoretical models [ 36, 37, 38] of the flux noise produced by weakly interacting electron 
spins located on the defects of the structure, these models differ in terms of the spin precise 
location. One possibility is that the spins are trapped at the interface between the supercon- 
ducting films and the insulating substrate 1381 . There is evidence 1 24, 26, 3] that some of 
the common insulating materials used for qubit fabrication have defects which affect the 
coherence times of superconducting qubits. The most likely candidates in our qubits are the 
substrate and the insulation under the wiring layer. We briefly review a number of tests con- 
ducted on our qubit to eliminate many other possible sources of decoherence. This leaves us 
with the conclusion that such TLFs are strong candidates for the decoherence seen here as 
well. 

One extraneous possibility is that the extensive filtering of external leads in not ade- 
quate. To check this, superconducting aluminum shunts were added to the flux bias coils 
allowing the bias flux to be trapped in a superconducting loop that included the shunt. This 
reduced the spectral density of any noise coupled to the qubit by the bias lines, as well as any 
damping due to the lines by two orders of magnitude but had no effect of the decoherence 
in the qubit !39l . The magnetometer was maintained at zero voltage until the qubit readout. 
By varying this zero voltage bias point the small coupling between the magnetometer bias 
lead and the qubit could be varied by several orders of magnitude and adjusted to zero (at 
least to first order). Again, no effect on coherence was seen from this variation. Coherence 
was also not affected by the repetition rate of the readouts, eliminating the possibility that 
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quasiparticle were being generated in the qubit by the readouts and affecting subsequent 
measurements. 

Back action from the magnetometer is another potential source of decoherence. How- 
ever, the magnetometer is operated in the switching current branch which is "silent" until 
it is time to perform the measurement. When the measurement is performed a voltage is 
established across the junctions creating quasiparticles which destroy the state of the qubit. 
However, for the technique used to measure the lifetime and the Rabi oscillations, the mag- 
netometer is left in the superconducting state until well after the state of the qubit was fixed. 
After the measurement, the magnetometer's bias current is ramped back to zero and left 
there long enough that the quasiparticles disperse. The magnetometer is designed such that 
its bias current does not couple strongly to the qubit. The flux and current biases can be set 
at a point where the derivative of circulating current with respect to the bias current goes to 
zero. At this point, fluctuations in the bias current do not be couple to the qubit to first order. 
However, there is no difference between operating at the optimum bias point and away from 
the optimum bias point in the measurements of the lifetime of the excited state. This implies 
that the magnetometer is not the major source of decoherence. 

We have also tested the quality of the junctions and Nb film used to fabricate our qubit. 
Patel et al. 1111 have found that the junctions co-fabricated with our qubit devices have a 
subgap resistance of 1 GQ or greater at 400 mK. This level of damping would give coher- 
ence times orders of magnitude longer than the values measured in Sec. [4] Critical current 
fluctuations, having a 1/f spectrum, have the effect of modulating the barrier height of the 
rf SQUID potential resulting in fluctuations of the energy spacing between coupled states 
causing decoherence. Recent measurements by Pottorf et al. 1401 of the level of 1/f crit- 
ical current fluctuations in junctions co-fabricated with our qubit devices have found this 
noise spectrum density to be orders of magnitude less than the commonly accepted value. 
Within the model of Van Harlingen et al. 1411 . these critical current fluctuations are orders 
of magnitude too small to explain the level of decoherence measured in our sample. 

The Nb thin film used for the roughly 0.5 mm flux loop in the qubit could potentially 
have excess loss despite being in the superconducting state. The losses at high frequencies 
(10's of GHz) were tested by measuring the Qs of coplanar waveguide resonators 1 42 1. The 
process for fabricating these resonators is similar to the process used to fabricate the qubit 
samples. Despite discovering losses at these frequencies that seem to depend on the fabri- 
cation process, the measured Qs (w 5 x 10 5 ) were sufficiently high to rule out these losses 
in the Nb films as a significant source of decoherence. We should note, however, that these 
resonator measurement were at too high a power level to exclude losses due to easily sat- 
urated TLFs. The effect of the substrate on decoherence has been studied by fabricating 
and measuring qubits and resonators on substrates of supposedly higher quality than the 
20 Qcm silicon substrates commonly used for our qubits. These included sapphire and 15 
kQcm silicon wafers, which showed much wider MRT peaks than those seen in Fig.flOlfor 
qubits fabricated on a 20 Q cm silicon substrate. Further investigation is necessary to check 
if the increased width is caused by the substrate and not by some unforeseen effect on the 
fabrication process. 

The SiC>2 layer between the junction layer and wiring layer is also a potential source of 
the observed decoherence. Most superconducting qubits with longer decoherence times are 
made using A1/A10 X /A1 junctions [43 . 1 1 using a two angle shadow evaporation process that 
has a thin layer A10 x for insulation. The NIST group uses a process the includes a deposited 
insulation layer. When they switched from SiC>2 to SiN x |3) they observed a significant 
increase in the decay time for their Rabi oscillations. 
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6 Conclusion 

With careful design of the measurement setup and a suitable fabrication process, it is pos- 
sible to observe coherent oscillations between two energy levels in an rf SQUID. Since the 
quantum coherence is extremely sensitive to external noise, extensive steps were taken to 
protect the rf SQUID. The sample was fabricated using a flexible, well characterized fabri- 
cation process capable of producing high quality Josephson junctions. The flux readout, a 
hysteretic dc SQUID magnetometer, was optimized to minimize the back action on the rf 
SQUID. Also the filtering and electronics were carefully designed to reduce the amount of 
external noise reaching the rf SQUID. Finally a setup for very fast pulsing of the bias flux 
and coupled microwaves was instituted to allow measurement of the coherent oscillations 
on a nanosecond time scale. 

The coherence times of the superposition between the ground and first excited state 
in the same flux state were measured and analyzed within the context of the known noise 
sources and calculations based on the parameters of the rf SQUID. There exists low fre- 
quency flux noise (0.14 m<I>o) and a short lifetime for the excited state (20 ns) that are not 
consistent with known noise sources. The coherence times are currently too short to be use- 
ful for quantum computation. However, it is likely that the decoherence is caused by issues 
relating to materials or fabrication and could be improved through a careful comparative 
study. 
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